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ABSTRACT 

A large number of spectroscopic studies have provided evidence of the presence of 
multiple populations in globular clusters by revealing patterns in the stellar chemical 
abundances. This paper is aimed at studying the origin of these abundance patterns. 
We explore a model in which second generation (SG) stars form out of a mix of 
pristine gas and ejecta of the first generation of asymptotic giant branch stars. We 
first study the constraints imposed by the spectroscopic data of SG stars in globular 
clusters on the chemical properties of the asymptotic and super asymptotic giant 
branch ejecta. With a simple one-zone chemical model, we then explore the formation 
of the SG population abundance patterns focussing our attention on the Na-O, Al-Mg 
anticorrelations and on the helium distribution function. 

We carry out a survey of models and explore the dependence of the final SG chem- 
ical properties on the key parameters affecting the gas dynamics and the SG formation 
process. Finally, we use our chemical evolution framework to build specific models for 
NGC 2808 and M4, two Galactic globular clusters which show different patterns in the 
Na-0 and Mg-Al anticorrelation and have different helium distributions. We find that 
the amount of pristine gas involved in the formation of SG stars is a key parameter to 
fit the observed O-Na and Mg-Al patterns. The helium distribution function for these 
models is in general good agreement with the observed one. Our models, by shedding 
light on the role of different parameters and their interplay in determining the final SG 
chemical properties, illustrate the basic ingredients, constraints and problems encoun- 
tered in this self-enrichment scenario which must be addressed by more sophisticated 
chemical and hydrodynamic simulations. 

Key words: stars: chemically peculiar - globular clusters: general - globular clusters: 
NGC 2808, NGC 6121 (M4) 



1 INTRODUCTION 

In the last decade, a large number of spectroscopic and pho- 
tometric studies have provided strong indication of the pres- 
ence of multiple stellar populations in globular clusters. 

The spectroscopic evidence comes from the observed 
spreads and anticorrelations between the abundances of light 
elements, not shown by the halo stars of similar metallici- 
ties. Specifically, all the clusters studied have shown an an- 
ticorrelation between Na and O and an anticorrelation be- 
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tween Mg and Al has been found in several massive clus- 
ters (cf . iGratton et al.l 12004 ) . The fact that these chemical 
anomalies have been found also in turnoff and subgiant stars 
(|Gratton et al.ll200ll ; lRarm'rez fc Cohenll2002r i indicates that 
they must be already present in the gas from which these 
stars formed. The antic orrelation between CN and CH bands 
observed in giant (e.g. iNorris et al.lll984T ) and turnoff stars 
indicates a difference in the nitrogen abundance and is con- 
sidered another fingerprint of the presence of multiple pop- 
ulations in globular clusters. 

An additional important constraint on the possible 
sources of gas from which second generation stars formed 
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is provided by the observed constancy of the total C+N+O. 
The only exception, so far, is NGC 1851 for wh ich a differ- 
ence of 0.6 dex is found among four red giants l|Yong et al,l 
[2009J) and for which the analysis of the subgiant branch 
suggests that the total C+N+O increases by about a fac- 
tor three betwee n the stars in the br ighter and those in the 
dimmer branch II Ventura et alj|2003 ). 

Photometric studies have also provided a number of im- 
portant results indicating the presence of significant differ- 
ences in the helium abundances of stars within individual 
clusters. Indeed, the presence of helium differences seems to 
be the only explanation for the multiple main sequences or 
spread in the main sequence (MS) observed in some clus- 
ters. These observations confirmed the prediction of the 
presence of a helium spread based on the complex mor- 
pholo gy of the horizontal branche s of some clusters made 
by |D'Antona fc Caloil |2004 120081 ). In extreme cases (the 
blue MS of w Cen an d of NGC 2808) a helium content 
Y=0.38-0.40 is in ferred (|Norrisll2004l ; lD'Antona et~ai1l2005l ; 
iPiotto et al]|2007l ). 

A recent spectroscopic survey bv lCarretta et al.l (|2p09c) 
has found that all the 15 clusters studied show the spectro- 
scopic evidence of the presence of multiple populations and 
that in all cases second generation stars represent a signifi- 
cant fraction (50-80 %) of a cluster population. This implies 
that matter forming more than half of each cluster stars has 
been processed through the hot CNO cycle and by other 
proton-capture reactions on light nuclei. 

In spite of the light elements variation, most GCs are 
mono-metallic objects, as far as abundances o f heavier el- 
ements are concerned (see iGratton et alj|2004l . for a recent 
review on this subject). Their heavy (Z >13) element metal- 
licity, usually represented by the ratio [Fe/H], is found to 
be e xtremely ho mogeneous from star to star in each clus- 
ter. ICarretta et al. (2009a) find that the upper limit to 
the scatter of iron is less than 0.05 dex in the 19 clusters 
they examine. This is consistent w ith previous deter mina- 
tions, e.g. for NGC 6752 fq =0.02 lYong et all I2005J ). and 
for the clusters examined by Kraft fc Ivans! ( 2003T ) (0.03 - 
0.10 dex). The width in color of the main sequence and/or 
of the red gian t branch also agree with the spectroscopic 
determination (|Suntzefll ll993T ). As the site of production 
of heavy elements (a-capture and Fe-group elements) are 
stars explo ding as core-collaps e or thermonuclear super- 
novae (e.g. I Wheeler et al.lll989r ). the star to star iron ho- 
mogeneity means that supernova ejecta do not affect the 
chemistry of the gas processed through the hot CNO cycle. 
On the other hand, significant iron and s-process elements 
spread is well known to be present in lj Cen (e.g., for a 
recent analysis, Ijohnson et al.ll2009l . and references therein) 
and more recentl y it has been confirmed in other massive 
clusters like M22 l|Marino et al.l | 2009l; |Pa Costa et ai1l200S , 
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1995; 



and references there in), M54 
Bcllazzini et al. 2008, and references therein), and Terzan 5 
( Ferraro et alj|2009l ). In smaller clusters like NGC 1851, the 
SG m ight be enriched in calciu m, according to Lee et al.l 
(2009), but see the dis cussion bv lCarretta et all l|2010h . 

In iD'Ercole et al.1 l|2008h (hereafter Paper I), we pre- 
sented a model for the formation and dynamical evolution 
of multiple populations in globular clusters (GCs). In partic- 
ular, we explored a model in which second generation (here- 
after SG) stars form out of the ejecta of the first generation 



(hereafter FG) stars. By means of hydrodynamic simula- 
tions we have shown that the ejecta of AGB (asymptotic 
giant branch) stars collect in a cooling flow into the clus- 
ter core, where they form a subsystem of SG stars initially 
strongly concentrated in the cluster innermost regions. By 
means of N-body simulations we have then explored the sub- 
sequent stellar dynamical evolution of the cluster focusing 
our attention on the early loss of FG stars, on the evolution 
of the ratio of the number of FG and SG stars and on the 
evolution of their relative spatial distribution and mixing. 

Paper I also included a preliminary discussion on 
the chemical abundances of FG AGB stars and on the 
possibility of reproducing the helium excess and other 
chemical anomalies observed in SG stars. Several the- 
oretical and observational papers have shown that di- 
lution of the FG ejecta with gas having the pris- 
tine composition is nee ded to explain the g e neral shape 
of th e anticorrelations (IPrantzos et al.l 120071; iBekki et al.l 



I2007J ; iD'Antona fc Ventural 120071 ; ICarretta et al.l l2009cT l 
If the SG anomalies are inputed to processing by Hot 
Bottom Burning (HBB) at the bottom of the con- 
vective envelope of massive AGBs, these models pro- 
duce a direct correlation between sodium and oxygen 
abundance in the processed matter. This direct correla- 
tion is unavoi dable, as found in all the relevant AGB 
computations (Denisscnkov & Hcrwig 2003; Hcrwig 2004; 
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and as explained in more detail in Section 2. While the to- 
tal sodium yield is dependent on many uncertain factors, 
it is certain that a small oxygen depletion (lower burn- 
ing temperature) goes together with a larger sodium abun- 
dance and a large oxygen depletion (higher burning tem- 
perature) is accompanied by a smaller sodium abundance. 
Any attempt to reverse this direct correlation, in the lim- 
its p rovided by available cross se ctions, has been unfruitful 
(e.g. IVentura fc D'Antonall2008br ). An obser vational hint in 
favour of the above arguments is given by ICarretta et al.l 
l|2009a ) who find a direct correlation between the minimum 
O and the maximum Na abundances of the 15 clusters they 
studied; such a result could not be understood if O and Na 
of the AGB ejecta were anticorrelated. 

From all the above, it is apparent that, in order to ex- 
plain the GC chemical patterns, it is necessary to introduce 
a form of dilution of the AGB ejecta, either with pristine gas, 
or with gas not showing the peculiar abundance patterns of 
the hot-CNO processed matter. 

The source of this pristine matter stil l requires a de- 
tailed understanding. ID'Ercole et al.l (J2008) proposed that 
initial asymmetries in the gas distribution allow to vent out 
the Type II supernovae (SNe II) ejecta along preferential di- 
rections, creating an "hour- glass" cavity and leaving some 
pristine gas in a torus at the outskirts of the cluster; at the 
end of the SN II activity this torus eventually collapses back 
into the cluster. In this scenario, the torus is not affected by 
contamination from supernova ejecta, and the iron content 
of the diluting matter is still the pristine one, consistently 
with the absence of iron differences between the FG and SG 
summarized above. While the evolution of the larger clusters 
is probably also affected by more prolonge d star formation 
and supernovae type la contamination (e.g. iMarcolini et al.l 
[2007J), the presence of slight metal increase in smaller clus- 
ters might be a possible indication that the "pristine" gas 
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in the falling back torus has been partially contaminated by 
the SN II ejecta. 

This paper is aimed at significantly expanding the ini- 
tial study of SG stars chemical abundances presented in Pa- 
per I and at exploring in detail the SG chemical anomalies 
emerging from a model in which SG stars form from a mix 
of AGB ejecta and gas with pristine chemical composition. 

We will refer to the scenario presented in Paper I. Con- 
sequently, we need to associate to the temporal evolution of 
the FG clusters the chemical composition of the ejecta of 
super-AGB and then massive AGB stars (hereinafter sim- 
ply referred to as AGB polluters) that successively provide 
the matter used for building up second generation stars. The 
chemical composition of the H-rich matter in the envelopes 
of these stars is affected by the hot CNO cycle processing at 
the basis of the convective envelope of these giants (hot bot- 
tom burning, hereafter HBB) and possibly also by the third 
dredge up (3DU) occurring after each thermal pulse. Stellar 
winds and planetary nebula ejection deposit the matter of 
these envelopes into the interstellar medium (ISM), where 
it can contribute to the formation of the second generation. 

Despite the significant effort devoted to the determina- 
tion of the yields of the massive AGBs and the numerous 
attempts to quantify the unc er tainties in these calculation s 
l| Ventura fc D'AntonalkoOSal lbl; iKarakas fe Lattanzidl2007l) , 
there are still large differences in the results obtained by dif- 
ferent groups and even by the same groups at different times 
and no general consensus has been reached. 

Our aim is to focus on the chemical properties of the 
stellar envelopes and explore the constraints on these prop- 
erties imposed by the spectroscopic data of SG stars in glob- 
ular clusters. We have addressed this problem by means of 
a simple one-zone model in which the ISM is supplied by 
the FG AGB ejecta and by the possible accretion of pristine 
gas, and is depleted by the SG star formation (hereafter 
SF). This simplified approach allows to easily control the 
evolution of the amount of the AGB ejecta and pristine gas 
involved in the SG formation process and to fully explore 
the chemical properties of the ISM and of the SG stars. 

The structure of the paper is the following. We start in 
sect. [2] with an overview of the chemical abundances of the 
massive AGB ejecta. We then illustrate our general assump- 
tions in sect. [3] and the model details in sect. [4] We then 
analyze the model in sect. [SI and compare our results with 
the data of NGC 2808 and M4 in sect. [6T] and sect. RH1 
respectively. Finally, we summarize our conclusions in sect. 

m 



2 EDUCATED GUESSES ON THE YIELDS 
FROM SUPER-AGB AND MASSIVE AGB 
STARS 

We start by re-examining all the uncertainties connected 
with the computation of the chemical yields of the AGB 
polluters, identifying the range of possible models produc- 
ing chemical patterns and trends consistent with those found 
in spectroscopic observational studies. We will adopt yields 
for m assive AGBs listed in Table 2 in I Ventura fc D'Antonal 
(2009) for metallicity Z=10" 3 (see Fig. []}, and in Table 
[l] modified and extrapolated as shown in the columns la- 
belled "n". For carbon and nitrogen, the abundances of 
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Figure 1. Time evolution of the chemical abundance of the ejecta 
of the dying stars adopted in th is paper (thick lines) and as given 
bv lVentura fc D'Antonal 1120091 ) (thin lines). On the top axis some 
masses of dying stars are reported in correspondence to their life- 
time. 



IVentura fc D'Antonal l)2009l ) are adopted without modifica- 
tion (for 6.5 and 9 Mq, the values of M=6.3 M© have been 
adopted). 

We list in the next subsections the problems relative 
to each individual chemical element, and how they are re- 
lated to the physical inputs (e.g. cross sections) or to the 
treatment of some important structural problems (e.g. the 
treatment of convection, mass loss, possible extra-mixing at 
the formal convective boundaries). We will mainly deal with 
the yields relative to the metallicity Z=10~ , [a/Fe]=0.4, 
for two reasons: I) this metallicity is adequate to describe 
most of the GCs for which data relevant for the study of 
the chemical anomalies are available; 2) this is the metallic- 
ity for which most theoreti cal mod els have been built (see 
IVentura fc D'Antona|[2008al lbl [2009l , for yields calculated for 
different values of Z) . 

2.1 The role of super— AGBs and their yields 

Most researchers focused on the problem of multiple stel- 
lar populations mainly when two clusters (w Cen and 
NGC 2808) were recognized to have a large fraction of main 
sequence stars (~25% in to Cen , ~15% in NGC 2808) in 
a "blue MS" which could only b e interpreted as a very he- 
lium rich MS jBedin et al.ll2004l; iNorris! 120041; IPiotto et all 



120051 ; iD'Antona et all 1200a IPiotto et all 120071 ) PI In both 
clusters, this extreme-helium population must also be very 
homogeneous, both in helium and metal content, otherwise 
the blue MS would not be detached from the rest of MS stars. 
Two other very massive clusters, NGC 6441 and NGC 6388, 
show an horizontal branch (HB) extended to very high-T e ff , 



1 Actually, IPiotto et all ||2007|) show that in NGC 2808 three 
different parts of the MS can be identified. 
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Table 1. Averaged abundances in the ejecta of massive AGB and super-AGB stars. 
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"Total evolutionary time until the AGB phase 
'Core mass at the b eginning of the AGB phase. 
c From ISiesd J2007ah. 
d From lVentura fc D'Antonal J2009I) . 
e Values adopted in this paper. 



and their RR Lyr variables have very long periods, indi- 
cating very high RR Lyr luminosity, a p eculiar occurrenc e 
at the high metallicity of these clusters (|Rich et al.lll9971) . 
Also in these cases, a very high helium population (including 
~ 10% and ~ 20% of stars respectively) provides the only 
reliable explanation for these features jCaloi fc D'Antonal 
120071 ; iBusso et al.ll2007l ; lYoon et al.ll200Sf ). 

In all cases, modelling the MS colours, or the HB mor- 
phology and the RR Lyr periods, the helium content of the 
most extreme population results to be Y~0.38-0.42. Values 
so high are extremely peculiar, particularly because the he- 
lium enrichment is not accompanied by metal enrichment (a 
value ~ 70 is estimated for the parameter AY/AZ). Mod- 
els of massive AGBs can reach Y values of ~ 0.35 in their 
envelopes, and therefore models based on the formation of 
this high helium population out of massive AGBs ejecta can 
not reproduce the h igh Y values suggested by observations 
IjKarakas et al.l l2006) . Similarly, models based on yields com- 
ing from either Type II or Type la supernovae (SNe la) are 
not successful (Romano et al] 120071: JMarcolini et all 120071; 



IChoi fc Yil [20081 ). IPumo et al] ([20081 ) noticed that a very 
high helium content (Y up to ~0.38) can be achieved al- 
ready just after the second dredge-up in the envelopes of 
the m ost massive s uper-AGB' stars, as shown by the mod- 
els by Sicss (2007a). Super- AGBs are the stars that reach a 
carbon oxygen core mass of ~ 1.05 Mq, and ignite carbon 
in conditions of semi-degeneracy. The final outcome of this 
burning is a degenerate oxygen-neon core, that can now tra- 
verse a "normal" thermal pulse (TP) phase, like the stars 
with C-0 core. A successful preliminary chemo-dynamical 
model, able to explain the helium distribution of the triple 
main sequence of NGC 2808, was then proposed in Paper I. 
We assumed that a "first" SG formed out of the ejecta of 
super-AGBs, so that this population has a very high Y, sim- 
ilar for all stars, while, as the SG formation phase proceeds, 
a "second" SG forms out of gas from massive AGBs mixed 
with pristine gas falling back onto the central regions and 
has a smaller helium content. This model produces a gap 
in the helium distribution of SG stars and can successfully 
reproduce the distinct populations identified in the color- 
magnitude diagram of NGC 2808. 

To make further progress in understanding the origin 
of SG stars and to test the viability of different hypotheses 



on the origin of the gas out of which SG stars formed, it is 
essential to explore whether all the chemical anomalies that 
characterize SG stars can be successfully reproduced. Unfor- 
tunately, only a few models of super-AGBs are available in 
the literature, and no low metallicity mass has been evolved 
through the whole TP phase. For the case of helium, we can 
rely on the helium abundance reached i n the envelope after 
the second dredge up, as computed by ISiesa (|2007al) . but 
we do not know how precisely the other abundances will be 
modified by HBB during the rest of the stellar life, until the 
whole envelope is lost to the interstellar medium. In addi- 
tion, even when the first models through the super-AGB TP 
phase will be available, their input physics must be carefully 
explored, as done for the massive AGB stars. Until then, we 
must make guesses on their yields. In this paper, we sim- 
ply propose that the yields from super-AGBs for O, Na, 
Mg, Al, are compatible with the patterns observed in the 
most extremely anomalous GC populations. For instance, 
the sum of the CNO abundances can remain close to the 
initial value during the thermal pulse phase of super-AGBs, 
as they have a limited efficiency of third dredge up, because 
of the weak helium luminosity reached during the thermal 
pulses (|Siessll2007bl ). We then extrapolate the yields com- 
puted for the lower AGB masses to the super-AGB regime. 
The extrapolation must satisfy two constrains: 1) super- 
AGBs have oxygen-neon core masses larger than the largest 
carbon-oxygen core masses of the most massive AGBs that 
do not ignite carbon. Consequently, we expect more extreme 
physical conditions in their convective envelopes, than in 
those of the most massive AGBs we have computed, and a 
more advanced nucleosynthesis due to proton-capture ele- 
ments (this assumption also depends on the effect of mass 
loss, and on the total duration of the thermal pulse phase); 
2) in addition, the results must not be in deep contrast with 
the observed abundance patterns. As for the effect of the 
third dredge up, in this work we will assume that dredge up 
in super-AGBs, if any, is less efficient than in less massive 
AGBs, as suggested bv lSiessI (|2007bl) . 
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2.2 Oxygen burning 

One of the most challenging issues in modelling the abun- 
dance patterns of SG stars in globular clusters is the SG low 
oxygen abundance. 

Population II stars are in general a-enriched, and thus 
oxygen-enri ched, due to the role of SNI I in the oxygen pro- 
duction fe.g. lMatteucci fc Greggiolll986r ); in order to deplete 
the abundance of oxygen, the gas from which SG stars are 
formed must have gone through the full CNO cycle (the ON 
chain) that occurs at temperatures of ~ 40 MK in the stel- 
lar interior, and needs temperature of Tb ce ~ 60 — 80 MK 
at the bottom of AGB convective envelopes (Tbce), both 
due to the much lower stellar densities of these regions, and 
to the shorter evolu tionary times of this phase (see, e.g. 
IPrantzos et al.ll200it l. 

Efficient ON burning occurs due to HBB in massive 
AGBs — and consequently in super-AGBs, due to the higher 
tempera tures reached at the bottom of their convective en- 
velopes. IVentura fe D'Antonal (|2005aT ) showed that more ef- 
ficient ON processing is favoured by more efficient convec- 
tion. However, the most massive AGBs with values of the 
metallicity (Z~ 10 -3 ) typical of the clusters showing very 
strong chemical anomalies, can not deplete oxygen by more 
than ~ 0.8 dex. Starting form a typical [O/Fe]=+0.4, the 

minimum oxygen abundance is [O/Fe] ~ -0.35 0.40. In 

most GCs the O-depletion is limited t o ,$ 0.8 dex, but we do 
know that gia nts in clusters like M 13 (|Sneden et al.ll2004l ) or 
NGC 2808 (|Carretta et alj|2008h show extreme stars with 
oxygen depletion down to a factor larger than 20. 

There are two options to explain the oxygen abundance 
of these extreme stars: 1) either the strong O-depletion is 
due to the higher Tbcc reached by super-AGBs; 2) or it 
is a consequence of in-situ mixing in red giants (these ex- 
treme anomalies have been so far found only along the giant 
branch) suffering deep extra-mixing due to their very high 
initial helium content. Indeed, in these stars, the molecular 
weight barrier that could forbid deep mixing at luminosities 
below the red giant bump is much sma ller than in stars with 
normal Y JD'Antona fc Ventur j|200^ . 

In both cases, if we wish to have an educated guess 
about the oxygen yield of super-AGBs, the choice is manda- 
tory: it must be smaller than for the maximum computed 
AGB mass, so we adopt [0/Fe]=-l.l for the progenitor 
largest mass 9 Mq. 



2.3 Sodium production 

Understanding the sodium abundance in the SG stars 
presents similar complexities. The sodium yield is provided 
by three different stages of the life of a massive AGB star. 
At each of these stages, different problems must be attacked 
to determine the yields, and may have different solutions. 
In addition, during these stages an AGB star loses mass, 
so the total yield will also depend on the mass loss rate 
in each of these phases. We have discussed these problems 
in a s eries of papers dedicated to the nucleosynthesis in 
AGBs (jVentura fc D'Antonall2008al . and references therein) . 
Here below we outline the phases during which the sodium 
changes: 

(1) The first phase during which sodium increases is when 
the second dredge up, before the star climbs up the AGB, 



mixes to the surface the 23 Na and 22 Ne produced inside the 
star during the previous H-burning phases. 

(2) Immediately after, the reaction 22 Ne(p,7) 23 Na takes 
place, due to the HBB, and 23 Na begins to be cycled back 
through the reaction 23 Na(p,a) 20 Ne. The most important of 
these two paths in order to determine the total 23 Na yield 
is the first one, 22 Ne(p,7) 23 Na, that, unfortunately, is un- 
certain by a factor ~2000 in the range of te mperature of 
interest for oxygen depletion (|Hale et alj|2002h . It is impor- 
tant to point out, however, that, once all the 22 Ne is burned 
to 23 Na via HBB, the sodium production does not increase 
even if we further increase the cross section of this reaction 
Ijlzzard et al.ll2007r i. We can keep mor e sodium only by low- 
ering the reaction rate 23 Na(p,a) 20 Ne (|Ventura fc D'Antonal 
120061 ). 

(3) The third stage by which sodium in the envelope 
changes is the third dredge up phase (3DU) following each 
thermal pulse. The smaller the initial mass, the more ef- 
fective the 3DU becomes. At each episode, the 3DU in 
fact mixes down into the helium layers the nitrogen ob- 
tained in the H-envelope due to HBB-CN cycling of the 
12 produced during the thermal pulse and partially dredged 
up. Consequently, primary 22 Ne, produced by the chain 
14 N(a,7) 18 F(/?+,i/) 18 0(a, 7 ) 22 Ne, is also dredged up from 
the helium inter-shell, and is converted into primary sodium 
by HBB. If the stars go through many episodes of 3DU, this 
third phase can increase the sodium abundance even by or- 
ders of magnitude, but it also leads to a strong increase in 
the total C+N+O yield (see sect. 12. 4[) which has not been 
observed in globular clusters. 



All the computations of AGB evolution in the liter- 
ature clearly show that, as the AGB mass increases, the 
sodium yield decreases; this is a consequence of the increas- 
ing Tbcc- A larger Tb co depletes more oxygen, and destroys 
sodium more easily by the 23 Na(p,a) 20 Ne cycle. In addition, 
in smaller masses the 3DU will dredge up more the 22 Ne 
and convert it into sodium. This implies that the massive 
AGB models, and probably the more massive super-AGBs 
as well, provide a direct O-Na correlation, and not the ob- 
served anti-correlation. We address this problem in sect. 15. ll 
and show the possible role played by dilution with pristine 
gas in turning this correlation into the observed Na-O anti- 
correlation. 

In order to reproduce the observed abundances of Na in 
the SG, we need to have the largest possible sodium yields for 
the masses in which the 3DU does not play a strong role, that 
are, as we will see from the C+N+O examination, M ^ 5 
Mq. Starting f rom the reference values of t he yields from the 
calculation of IVentura fc D'Antonal (J2009J) (see their Table 
2 and Table 1), there are two key ingredients which are still 
affected by significant uncertainties and could lead to larger 
sodium abundance in the massive AGB envelopes: 



(i) the 20 Ne(p,7) reaction rate: a 50 % increase in 
the 2 Ne(p,7) rate is well within the NACRE (|Angulo et all 
Il999l ) reaction rate uncertainties and, keeping all other pa- 
rameters fixed, will lead to a 0.15 dex increase in the sodium 
yield, because the equilibrium of the Ne-Na chain is shifted 
towards 23 Na; 
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(ii) the 20 Ne abundance: the initial 20 Ne abundance is 
not well knowrQ and it may also vary from cluster to clus- 
ter (sect. [6}. In our reference computation, the Ne abun- 
dance is already a-enhanced, as for oxygen and the other 
a-elements, by 0.4 dex with respect to the solar ratios. If 
we increase this initial 20 Ne by another factor 2, the final 
[Na/Fe] is increased by ~0.2 dex (for the important masses 
from 6.3 to 5 Mq). Owing to our numerical tests, in a first 
approximation we can assume: 



S[Na/Fe] ~ S[ 20 Ne/Fe] 



(1) 



If we adopt both the enhanced reaction rate and the 
larger 20 Ne abundance, the combined effect of these two 
assumptions leads to [Na/Fe] = [Na/Fe] s tandard+0.35. In our 
calculations the increase of the sodium yield is calibrated 
by comparison with the observational data by +0.3dex (see 
Tableland Fig.[TJ for the masses up to 6 Mq . For the super- 
AGB stars, we decrease the sodium content to [Na/Fe]=0.45, 
based on the hypothesis that these stars will be more effi- 
cient in destroying sodium. We warn that this suggestion will 
need confirmation from real computation of super- AGB evo- 
lution. Once we have assumed a "standard" table of abun- 
dances with increased sodium, we still can use equation [T] 
to model clusters with abnormally larger a-elements abun- 
dances (see the case of M4 in sect. I6.2JI 



2.4 The total CNO abundance in the SG 

For C, N and O we will u se the yields from Table 2 of 
I Ventura fe D'Antonal i|2009f ) without any adjustment. For 
the masses from 6.5 to 9 Mq, the super-AGB regime, we 
adopt the C and N envelope abundances of the most mas- 
sive computed AGB model of 6.3 Mq. The O abundance in 
this mass range, instead, has been calibrated to fit the abun- 
dance of the O-poor stars present in NGC 2808, as discussed 
in the previous subsection (cf. Fig. [l|. 

As most of the SG stars in GCs do not show a sig- 
nificant increase of the total CNO |lvans et al.lll999l . feoOll ; 
ICarretta et al.ll2005l ). the matter out of which the SG forms 
can not come from stars in which the effect of the 3DU is 
relevant. A reasonable assumption based on the data is that 
the total CNO in most clusters so far examined is within 
a fac tor two of the CNO o f the F G. If we adopt the yields 
from I Ventura fe D'Antonal (J2009J) as our reference values, 
and assume that SG stars form only from pure AGB ejecta, 
the limiting mass for the SG formation is then ~ 5 Mq, and 
the limiting age ~100 Myr (time evolution of the 5 Mq, see 
Table [TJ. If, on the other hand, the matter from which SG 
stars form is diluted with pristine gas, as shown to be nec- 
essary in sect.O the limiting age for SG formation can be a 
bit larger. 

An interesting indication concerning possible variations 
in CNO comes from recent spectroscopic observations of a 
few red giants in NGC 1851. These observations have shown 



2 In recent years, an increase in the 20 Ne abundance for the 
Sun has been suggested, as it could produce the opacity increase 
necessary to obtain a better fit of the solar oscillation spectrum, 
once the iron and oxygen content are reduced following the results 
of the 3D analysis of the solar spectrum ( for a summary of this 
complex problem, see lAsplund et al.l l2009). 



that the C+N+O abundance of the stars obser ved varies 
by up to a factor 4 (0.6 dex, lYong et al.1 20091). In addi- 
tion, new HST photometric data i Milone et al.l 120081 ) have 
revealed a split subgiant branch. Two interpretations are 
possible: either t he stars in the lowe r subgiant branch are 
older by ~ lGyr (JMilone et al.ll2008h . or they are SG star s 
with a larger total C+N+O content (jCassisi et al.1 120081 ). 
IVentura etahl (|2009| N ) show that the total C+N+O must be 
about a factor three larger in the lower subgiant branch 
stars. They argue that the matter from which the SG formed 
should have been diluted at 50% with pristine matter hav- 
ing approximately the composition of the ejecta of ~ 4 Mq 
AGBs. With this choice, they can explain at the same time 
the C+N+O abundance and the lack of a significant helium 
enrichment, that would be reco gnized in the lumin osity of 
the blue horizontal branch stars (ISalaris et alj|200a ). In our 
scenario, the SG in this case should have formed from mat- 
ter collected in the clus ter until an age of ~ 150 Myr, at the 
evolution of the 4 M Q l|Ventura et alj|2009h . 

In the construction of possible models for the SG, 
then, the C, N and O composition of the SG is an 
important parameter, and it is also important to ob- 
tain further observations of its range of variability for as 
many clus ters as possible. We point out that the AGB 



models of IVentura fe D'Antonal (|2009T ) adopted here have 
small er CNO e nhancemen t tha n other published models 
(e.g. Karakas fe Lattanziol I2007T ). as they have a smaller 
number of TPs and 3DU episodes. This is due both to the 
mass loss rate adopted and to the use of an efficient con- 
vection model that provides a larger stellar luminosity (and 
consequently a larger mass loss rate) for a fixed core mass. 
Should we adopt models with a more efficient 3DU, or a 
larger number of TPs, the formation of the SG must occur 
at shorter age (larger evolving initial masses) both for the 
"normal" SG clusters (no CNO enrichment) and for clusters 
like NGC 1851, where the CNO enrichment is present. 

2.5 The Magnesium and Aluminium yields 

In some (massive) GCs only, also an anticorrelation between 
magnesium and aluminium is revealed. The depletion in to- 
tal Mg is restricted to 0.2 - 0.4 dex. Generally, the stars with 
the highest degree of oxygen depletion are also strongly en- 
riched in Al, up to the values [Al/Fe]~ 1 — 1.2 found in the 
red giants of M13 a nd M15 (for a sum mary of the data of 
the last decade, see llvans et al.lll999r ). The anticorrelation 
shows that the matter of the SG has been subject to the 
Mg-Al chain, that, favouring proton captures by the heavy 
isotopes of magnesium, eventually leads to the synthesis 
of 27 A 1 (|Denissenkov fe Herwij|2003l : IVentura fe D'Antonal 
l2005aT ). In models with many third dredge up episodes, pen- 
etration of the convective envelope into the layers touched 
by 3a burning may bring to the surface 25 Mg and 26 Mg 
synthesized via a captures on 22 Ne nuclei. These isotopes 
produce 27 A1 by proton capture, onc e the H-shell is reac- 
tivated l|Ventura fe D'Antonal l2008bT ). The recent analysis 
of Mg and Al abundance s from UVES spectra in 18 GCs, 
by ICarretta et al.1 (J2009bj), shows however large cluster to 
cluster variations in the Al-Mg anticorrelation, and that in 
most cases the Mg variation is very limited, even for stars 
with [Al/Fe]~0.8. They also show that the maximum Mg 
abundance and the minimum Al abundance (that should be 
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indicative of the composition of the FG) may be significantly 
different in different clusters -even by an order of magnitude 
for Al. Where the Mg abundance is larger than the typical 
[Mg/Fe]+0.4 assumed for the a-enhanced composition, this 
may reflect different overabundances in the matter forming 
the FG, and we must take it into account, both for magne- 
sium and for the other a ele ments (see sect. 16.21). 

In their theoretical study. I Ventura &: D'Antonal (|2008bl . 

2009) adopted the NACRE upper limits for the aluminium 
production from proton captures on 25 Mg and 26 Mg, in order 
to obtain very significant Al production in their massive 
AGB models. Their aim was to reproduce the extended Mg- 
Al anticorrelation present in the MI3 and MI5 red giants. 
Their results, however, are in contrast with the abundance 
patterns in M4 (NGC 6121) in which the maximum [Al/Fej 
is ^0.8, and the min imum is [Al/Fe]~ 0.6 divans et al.ll 19991 ; 
iMarino et a l. 2008a), so that this cluster does not show great 
Al variations, if any. 

As the large Al abundances of NGC 2808 seem to re- 
fer only to very oxygen poor stars, we can attribute them 
to the super-AGB evolution. Consequently, we decided to 
restore the above cross section to their standard values, 
and assume that significant Al production is limited to the 
super-AGB masses only. For the temperatures of interest 
here (T ~ 10 8 K), the NACRE recommended values for the 
25 Mg(p,7) 26 Al and 26 Mg(p,7) 27 Al reactions are smaller, re- 
spectively, by a factor of 2 a nd 4 than the rates adopted by 
I Ventura fc D'Antonal (|2009r l . The aluminium production is 
consequently reduced: its average content in the ejecta of 
massive AGBs of 5 and 6Mq with metallicity Z=0.001 di- 
minishes to [Al/Fe]=0.7 and 0.62 when the recommended 
cross sections for the proton capture reactions by magne- 
sium isotopes are adopted, to be compared to the published 
values of [Al/Fe] = 1.02 and 1.04. The modifications to the 
Mg and Al abundances are provided in Table Q] 



3 GENERAL ASSUMPTIONS OF THE MODEL 

The general assumptions driving the present models are the 
same adopted in Paper I, and we recall them here briefly. We 
assume that the FG stars are already in place and have the 
same chemical abundances of the pristine gas from which 
they form. The SG stars form from AGB ejecta. In order to 
have a current similar number of SG and FG stars, the FG 
population must have been initially about ten times more 
massive so as to supply enough mass in AGB ejecta to form 
a sizable SG population. This conclusion remains substan- 
tially valid even if, during its evolution, the GC accretes 
some e xternal pristine gas. A very large initial FG is invoked 
also by iRenzinil (|2008l ) on the basis of the large amount of 
helium observed in some multiple populations of very mas- 
sive GCs (see sect. 2.1). 

We assume that only SG stars with M ^ 9 M© can 
form (cf. Paper I). More massive stars expl ode as SNe II, 
but th e very small spread of iron in GCs (cf. ICarretta et al.l 
2009a) indicates that the SG stars hardly have been polluted 
by SNe II (we discuss this point in the Introduction and in 
sect. [7}. Moreover, as shown in sect. [6] the AGB pollution 
takes tens of Myr to produce the observed chemical patterns; 
as the first SNe II starts to explode after ~ 5 Myr, they 
would vent away all the gas, preventing the formation of a 



substantial amount of SG stars Q Given a Kroupa initial 
mass function (IMF) with an upper mass of 9 M© , the mass 
of the long-lived (M < 0.8 M ) SG stars is 54% of the total 
forming SG population. 



4 COMPUTATIONAL METHOD 

The computational framework adopted for this study is 
based on a one-zone model for globular clusters in which the 
evolution of AGB ejecta and pristine gas - from which SG 
stars form - is regulated by a few parameters. The variation 
of these parameters can be easily controlled and this allows 
us to explore their relative importance on the resulting SG 
chemical properties. 

The evolution of the gas and the process of SG forma- 
tion is modelled by the following equations which are inte- 
grated by means of a fourth-order Runge-Kutta method: 



... . ... pit) 

pit) = ap*,PG + Ppr(i) - v— — 

p\{t) = a k p*,FG + (3 k p pl (t) - v^ 



Ppr(*) = ~^exp- (i-iac)/T ~V— 



P*,SG(t) = V 



pW 

fsf 



(2) 
(3) 
(4) 
(5) 



Here below we discuss each of the above equations: 

Eq. rj| - This equation describes the time evolution of the 
total gas density p. The first term in the right hand side 
(RHS) represents the increment of p due to the gas lost by 
the dying stars. In principle, we would take into account the 
mass return from the FG stars as well as from the newly 
formed SG stars. These latter, however, are much less than 
the former, and we neglect their contribution. We thus con- 
sider only the density p*,fg, and the specific mass return 
rate is given by a = p*,fg /p*,fg; this rate is computed fol- 
lowing [Ciotti_et_alJ l|199ll ) (see their eq. 3), adapted to the 
Kroupa IMF adopted here. The second term in the RHS of 
eq. [5] takes into account the gas increment due to the ac- 
cretion of the pristine gas whose density is p pr . Finally, the 
last term describes the reduction of gas because of the star 
formation, and is characterized by the SF efficiency v and 
the SF timescale tst. The former is a free parameter sum- 
marizing our uncertainties about star formation and with 
values in the range < v < 1; the latter coincides with the 
free-fall time, and is thus proportional to (p~ 2 jj. 



3 On the contrary, massive stars in the FG must be present for 
two reasons. First, as the FG is initially nearly ten times larger 
than today (in order to deliver enough ejecta to form the observed 
amount of SG), SN explosions allow the cluster to get rid of the 
mass excess in a time shorter than the Hubble time (cf. Paper 
I); second, Extreme, very O-poor stars present in some cluster 
(e.g. NGC 2808) can arise only from a pure AGB ejecta, and this 
could not be possible without the previous action of the FG SNe 
II clearing initially the GC of its residual pristine gas. 

4 In a full hydrodynamic treatment, as in Paper I, £ s f is defined 
as the maximum between the local cooling time and the local 
free-fall time, but the first one is usually much shorter 
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Eq. - Besides the total gas density, it is necessary to 
evaluate the density of the single chemical elements we are 
interested in. In our models we take into account the abun- 
dance evolution of nine elements: He, Li, C, N, O, Na, Mg, 
Al, Fe. Equation [3] represents a system of kmax = 9 equa- 
tions, each modelling the evolution of the density p^ h of the 
k — th element. In these equations a(t) and /3 k represent 
the mass return rate of the k — th element in the FG, and 
the mass fraction of the k — th element in the pristine gas, 
respectively. The latter is determined by the assumed initial 
chemical composition of the pristine gas, while the values of 
a(t) are obtained at each time by the interpolation of the 
values given by our adopted ejecta abundances (cf. Fig. [l] 
and Table [1| . 

Eq. p] - This equation describes the evolution of the pris- 
tine gas density p pr . The two terms on the RHS model the 
time evolution of the pristine gas accretion (first term) and 
the SG star formation (second term). The gaussian form 
adopted to model the pristine gas accretion is arbitrary, but 
allows us to easily explore different accretion scenarios. The 
time £ ac regulates the time at which the maximum accre- 
tion rate occurs, and the time r controls how the accretion 
evolves: a small value of r corresponds to a sudden accre- 
tion, while a large value gives rise to a more gradual growth. 
Finally, the constant po.pr controls the amount of available 
pristine gas. The three constants t a c, t and po.pr must be 
regarded as free parameters of our models. 

Eq. \5\- Finally, equation [5] regulates the growth of the SG 
stellar population. 

In addition to po, P r, f, iac, t, one more parameter must 
still be added: the time t en d at which the simulations stop. It 
regulates not only the amount of SG stars formed, but also 
their overall chemical distribution, as the chemical abun- 
dances of the FG ejecta vary in time (cf. Fig. [1} . Its value is 
constrained by the observational evidence (discussed in sect. 
12. 4|) that the sum of CNO eleme nts varies no more than a 
factor two among the GC stars (I Cohen fe Melendez 2005) 
(or by a factor three in NGC 1851. see lVentura et aiT l2009): 
this implies that the SG progenitors must be AGB stars 
with masses larger than ~5 Mq (~ 4 Mq in the latter case) 
and that the SG formation can not continue longer than the 
lifetime of these stars, that is i cn d < 100 Myr (or < 150 in 
the latter case) Q (cf. Fig. Q] for the behaviour of C+N+O 
with time). 

Finally, we introduce the ratio x = pS,sG/P*,tot be- 
tween the SG stars and the total nowadays alive stars 
p*,tot = p" fg + P* sg- ^s current value depends not only on 
the duration and efficiency of the SF of the SG stars, but 
also on the subsequent loss of FG stars (see Paper I). We 
point out that this parameter does not enter directly into 
the models (it is in fact absent in eqs. I2I5|) ; its value is in- 
ferred a posteriori for a more realistic fit of the data. We will 
return on this point in sect.EJ when comparing our models 
to specific GCs. 

In the end, we summarize in Table[2]the key parameters 
(besides the chemistry of the AGB ejecta) characterizing our 
models. 



5 A possible mechanism able to clear the ISM in a globular clus- 
ters in this lapse of time is given by the explosions of type la 
supernovae (see Paper I). 



Table 2. Parameters characterizing the models 

t ac time of maximum accretion of pristine gas 

t time scale of the pristine gas accretion 

t cnt j time marking the end of the SG formation phase 

P* FG density of the FG stars 

Po.pr density regulating the amount of available pristine gas 

v SF efficiency 

x the ratio between the nowadays alive SG and total stars 
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Figure 2. Top-left panel: chemical evolution of the RM in the 
plane [Mg/Fe]-[A1/Fe]. Bottom-left panel: chemical evolution of 
the RM in the plane [0/Fe]-[Na/Fe] plane. In these panels the 
dot-dashed line represents the ISM chemical pattern of the model 
without accretion of pristine gas (RMNOPG), while the solid line 
holds for the RM, in which the accretion occurs; the open circles 
represent a statistical sampling of the SG stars. Top-right panel: 
helium distribution function of the SG stars. Bottom-right panel: 
evolution of the total amount of ISM (bold solid line). AGB ejecta 
(solid line), pristine gas (long-dashed line) and SG stars (short- 
dashed line). 



5 MODEL ANALYSIS 

5.1 The reference model 

To illustrate how our model works, we focus 
our attention on a model with the following pa- 
rameter set (tac,7,tend,7,T7,pO,pr, p*,FG,t/,x) = 
(6.5,10,0.6,15,150,0.5,1), where the times are expressed 
in 10 7 yr and the densities in M© pc -3 . Note that we 
here are interested only in the evolution of the SG stars, 
and we set x = 1. We checked that the current choice of 
parameters gives (in absence of the accretion of pristine 
gas) the same temporal profile of the gas and of the SG 
stars of the reference model of Paper I (cf. its Fig. 2); for 
this reason, we'll refer to the one zone model described in 
this section as the reference model (RM). We consider an 
amount of available pristine gas one tenth of the FG mass, 
nearly a factor two larger than that of the ejecta delivered 
after ~ 10 8 yr. 
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Figure 3. Time evolution of the chemical abundances in the ISM 
of the RM. 



The main properties of this model are shown in Fig. [2] 
The two left panels illustrate the chemical path of the model 
in the planes [Al/Fe]-[Mg/Fe] and [0/Fe]-[Na/Fe]. The top- 
right panel provides the helium distribution function (here- 
after HDF) of the SG stars, while the bottom-right panel 
shows the time evolution of the amount of several quantities. 
Initially, only gas returned by the FG stars is present (thin 
solid line) . Later on, the total gas content (thick) first expe- 
riences a "jump" as a consequence of the pristine gas (long- 
dashed line) accretion, and successively decreases because of 
the star formation. The amount of SG stars (short-dashed 
line) shows, as expected, a rapid increase in conjunction with 
the accretion. 

Figure El illustrates the abundance evolution of different 
elements in the ISM of the RM. A comparison of Fig. [T] and 
Fig. El clearly shows that, as anticipated in sect. 12.21 the 
presence of pristine gas is necessary in order to reproduce 
the O-Na anticorrelation observed in all the GCs. In Fig. 
[T]the O and Na abundances in the FG stellar ejecta grow 
simultaneously, and no anticorrelation between these two 
elements can be established; on the other hand, Fig.[3]clearly 
indicates how the pristine gas accretion produces opposite 
slopes in the time evolution of the O and Na abundances in 
the ISM, thus producing an anticorrelation similar to that 
observed in globular clusters. 

This point is illustrated even more vividly in the 
bottom- left panel of Fig. [2] This panel includes the chemical 
path of the gas in absence of any accretion (i.e. p pr — 0), 
represented by the dot-dashed line. This model (hereafter 
the RMNOPG model) starts from the O-poor end of this 
line since the ejecta of the most massive AGBs are assumed 
to be particularly deprived of oxygen (see Fig. Q] and the 
discussion in sect. [2]). As gas from less massive AGB stars 
contributes to the ISM both the abundances of Na and O 
increase giving rise to a correlation rather than the observed 
anticorrelation. The path followed by the RM when the ef- 
fect of pristine gas accretion is included clearly illustrates 



the role of this gas in establishing the Na-O anticorrelation. 
Since the ISM is initially composed only by AGB ejecta, the 
RM track initially coincides with that of the model with no 
pristine gas. At t ~ 55 Myr, when [0/Fe]> —0.5, the amount 
of the infalling pristine gas becomes substantial, and the RM 
moves "downward" at increasing values of [O/Fe] , giving rise 
to a [0/Fe]-[Na/Fe] anticorrelation. This motion is not uni- 
form because the pristine gas accretion suddenly increases 
(decreases) the abundance of oxygen (sodium), "pushing" 
the model in short time towards higher (lower) values of 
[O/Fe] ([Na/Fe]) (cf. Fig. [3]). As the amount of pristine gas 
in the ISM increases, the trajectory moves toward its mini- 
mum located at [O/Fe]=0.37 and [Na/Fe]=0.1 dex, "above" 
the point identified by the ratios [O/Fe]=0.4 and [Na/Fe]=0 
of the pristine gas. Successively, the accretion rate of the 
pristine gas reduces becoming lower than the SFR, and the 
relative amount of this gas in the ISM decreases more and 
more. The ISM composition evolves back toward the chem- 
ical abundance of the AGB ejecta and the chemical path of 
this model tends again toward that of the RMNOPG model. 

It is important to point out that the ISM trajectory in 
the Na-0 plane is not sufficient to predict the distribution 
of SG stars in this plane. In fact, such a path is, in a sense, 
similar to a stellar Color-Magnitude diagram in which some 
regions are more populated than others owing to the longer 
time spent there by the stars. In our model the amount 
of SG stars along a certain portion of the chemical trajec- 
tory depends on the rate of evolution of the ISM chemical 
abundance as well as on the SFR (which is proportional to 
p 3/2 ). In order to show how the [O/Fe]- [Na/Fe] and [Mg/Fe]- 
[Al/Fe] diagrams are actually populated by the SG stars, in 
Fig. [2] we have plotted twenty open dots representing a sta- 
tistical sampling of the SG stars forming along the chemical 
path. Consider, for instance, the diagram [O/Fe]- [Na/Fe]: 
only a few stars populate the two branches of the curve, 
while most of them form close to its cusp; the cusp coin- 
cides with the maximum amount of ISM and the consequent 
maximum SFR. Although the time spent by the model along 
the two branches is similar (cf. Fig. El and the bottom-right 
panel of Fig. [2j , the descending one is more populated by 
stars because the ISM density and the SFR are larger at this 
stage. 

Finally, the top-right panel of Fig. [2] shows the HDF. 
A small number of stars (10%) have an helium abundance 
Y > 0.35. These stars coincide with the most O-poor stars 
present in the [O/Fe]- [Na/Fe] diagram. The stars giving rise 
to the peak around the range 0.24 < Y < 0.29 form when 
the maximum amount of pristine gas is present in the GC, 
and correspond to the stars populating the region around the 
cusp. The tail at Y ~ 0.3 is composed by stars formed during 
two nearly chemically "symmetric" phases, just before and 
after the occurrence of the bulk of the pristine gas accretion 
(cf. Fig. [3])). 



5.2 Exploring the parameter space 

We now explore the dependence of the model on the chosen 
set of parameters through variations of the chemical paths 
in the [0/Fe]-[Na/Fe] plane and variations of the HDF. We 
consider several models, each differing from the RM in the 
value of a single parameters. All the models are listed in 
Table El 
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Table 3. Models adopted for the exploration of the parameter space. 



Models 


^ac,7 


TY 


^cnd,7 


Ppr 


P..FG 


V 




panel and style c 


RM 


6.5 


0.6 


10.0 


15 


150 


0.5 




all, solid, black 


RMNOPG 








0.0 








all, dot-dashed, green 


RMMPRL 








1.5 








a, dashed, blue 


RMMPRM 








45 








a, dotted, red 


RMTACL 


5.0 














b, dashed, blue 


RMTACM 


8.0 














ft, dotted, red 


RMTAUL 




0.3 












c, dashed, blue 


RMTAUM 
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5.2.1 The [0/Fe]-[Na/FeJ diagram 

The six panels in Fig. [4] illustrate how changes of each single 
parameter affect the RM behaviour in the [0/Fe]-[Na/Fe] 
diagram. The RM path (solid line) is replicated in every 
panel to facilitate the comparison with the different modi- 
fied models. For the same reason we also include in all the 
panels the dot-dashed line representing the behaviours of 
the RMNOPG (see Fig. |2J. We now discuss these panels. 

Panel a - In this panel the RM is compared with two sim- 
ilar models, but with p pl = 4.5 M Q pc" 3 (RMMPRL, blue 
dashed line) and p pl = 45 M Q pc" 3 (RMMPRM, red dot- 
ted line). The O-poor branch of the RMMPRL substantially 
overlaps the analogous RM branch as long as the amount 
of accreted pristine gas is relatively low in both models. 
The rising branch, instead, is shifted toward lower values of 
[O/Fe] . This shift is a consequence of the smaller amount of 
pristine gas, and its lower influence on the chemical abun- 
dance of the ISM; the dashed line tends to return more 
rapidly toward the dot-dashed one. The reduced value of p pr 
affects the position of the cusp, which occurs at [O/Fe] =0.3 
and [Na/Fe]=0.2 earlier than in the RM model. 
Also in this case the stars (open blue squares) form mainly 
around the cusp, but the two branches are more populated. 
This is due to the combined effect of the lower amplitude 
of the gaussian accretion rate and the dependence of the 
SFR on p 3 , the result is a less marked peak in the SFR 
occurring at the cusp. 

Similar considerations apply to the RMMPRM. While the 
descending branch overlaps again with those of the previ- 
ous models, the cusp of the red dotted curve is located at 
a value of [Na/Fe] nearly 0.2 dex lower than that of the 
blue curve, and the rising branch is shifted toward higher 
values of [O/Fe]. Finally, in this model the stars (open red 
triangles) are essentially all concentrated around the cusp. 

Panel b - The models represented here are similar to the 



RM, but with tao,7 = 5 (RMTACL, blue dashed line) and 
t ac ,7 = 8 (RMTACM, red dotted line). The O-poor descend- 
ing branch of the dashed line is essentially similar to that of 
the RM, but slightly shifted toward lower values of [Na/Fe]; 
this is a consequence of the earlier accretion of pristine gas. 
The rising branch is instead longer and strongly shifted to- 
ward lower values of [O/Fe] because at this stage the amount 
of pristine gas eroded by the SF is larger than in the RM. 
The dashed line tends therefore to reconnect earlier to the 
RMNOPG chemical pattern. 

Given the delay in the gas accretion, the RMTACM follows 
the RMNOPG track for a longer interval, then quickly con- 
verges toward the cusp. The dotted rising branch is shorter 
and shifted toward higher values of [O/Fe] for opposite rea- 
sons relative to the RMTACL. 

Because of the delay in the accretion, the RMTACM de- 
scending branch is more populated by stars than the rising 
branch, while the opposite happens for the RMTACL. Once 
more, this is due to the fact that the SFR grows with the 
gas density. Note the gap in the stellar distribution along 
the descending dotted branch (open triangles). Initially, the 
RMTACM has enough time to form a substantial amount of 
AGB ejecta before the accretion becomes significant. A non 
negligible fraction of O-poor stars can thus form. The sub- 
sequent accretion moves the model toward the cusp very 
rapidly, and very few stars can be produced during this 
stage. 

Panel c - In this panel the models with Tr — 0.3 (RM- 
TAUL, blue dashed line) and r 7 = 1.2 (RMTAUM, red dot- 
ted line) are compared. Smaller values of r correspond to a 
more rapid accretion. The two lines are rather similar, but 
with the dotted rising arm being shorter and shifted toward 
O rich values; in this case, in fact, the longer duration of 
the accretion produces a more persistent influence on the 
RMTAUM which thus shows some delay in its tendency to 
reconnect with the RMNOPG line. 
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reported in all panels. The different symbols represent statistical samples of the SG stars in different models. 



The distributions of the stars along the dashed and dotted 
lines are also similar. However, a prevalence of O-poor stars 
in the RMTAUL (open squares) is present; owing to the 
shorter extent of the accretion, a larger amount of AGB 
ejecta can be initially collected without being significantly 
diluted by the pristine gas. 

Panel d - The effects of the SF efficiency are investigated 
by models with v — 0.1 (RMEFFL, blue dashed line) and 
v = 1 (RMEFFM, red dotted line). As usual, the descend- 
ing branches essentially overlap. The rising RMEFFM arm, 
instead, is more extended and shifted toward lower O abun- 
dances because the larger SF efficiency leads to a more rapid 
consumption of the ISM in general, and of the pristine gas 
in particular; once the bulk of the accretion is finished, the 
fresh AGB ejecta is diluted by a lower amount of the pris- 
tine gas still present in the ISM, and the RMEFFM chemical 
pattern moves more quickly toward the RMNOPG track. 

Panel e - In all the previous models, we assumed £ en d,7 = 
10. Obviously, the duration of the SG formation phase may 
be different, and we studied the influence of such a dura- 
tion on the chemical path of the SG stars. In this panel we 
show a model with t cn d.7 = 7 (RMTNDL, blue dashed line) 
and a model with t cnd , 7 = 13 (RMTNDM, red dotted line). 



As expected, the chemical paths coincide, but with the dot- 
ted one obviously longer than the dashed one. As the two 
models are identical, the number of stars formed along the 
common tracks is the same, and one could expect the same 
distribution of triangles and squares on them. Instead, it is 
apparent that the triangles are less than the squares on the 
first branches; this is due to the fact that the number of 
symbols indicating the stars is not proportional to the total 
number of stars actually formed; rather, for every model we 
adopt twenty symbols to represent the total number of SG 
stars in it, no matter if this number is different in different 
models. In particular, the two models shown in this panel are 
identical, and the same number of stars is formed along the 
descending branches; however, the RMTNDM lasts longer, 
and many more stars are formed along the rising branch. 
For this reason more triangles are distributed on this branch 
and less on the descending one, where the squares are more 
numerous. 

Panel f - In this panel we show the behaviour of two mod- 
els having both the values of p*,pg and p pl ten times lower 
(RMFGDL, blue dashed line) and ten times higher (RM- 
FGDM, red dotted line) than in RM. In sect. [5T] we have 
justified our choice of p,,FG — 150 Mq pc -3 ; this choice, 
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however, although reasonable, remains arbitrary. We expect 
that our results depend on a scaling of p*,fg and p pr , even if 
their ratio, and thus the relative amount of dilution due to 
the pristine gas, remains the same; in fact, a lower (higher) 
density of p*,fg produces a lower (higher) density of the 
AGB ejecta which, together with a lower (higher) p pr , de- 
termine a lower (higher) SFR and a different evolution. 
This difference is apparent in the present panel. As usual, 
also the two models shown here have very similar descending 
branches, but the RMFGDL cusp occurs at a higher value 
of [Na/Fe]. This is because, owing to the lower SFR in this 
model, less ejecta are initially subtracted by the star forma- 
tion, and the dilution by the pristine gas is less effective. 
Successively, however, always because of the lower SFR, the 
pristine gas is depleted more slowly, remaining longer in the 
ISM and leading to a shorter rising branch. Opposite argu- 
ments hold for the RMFGDM, whose cusp is "deeper" and 
whose rising branch is longer. Note also that this arm is 
strongly bent "to the left" because, given the larger amount 
of ejecta (relative to the pristine gas) in this case, this model 
tends more quickly toward the RMNOPG. 



5.2.2 The helium distribution function 

Another crucial output of the models are the stellar HDFs, 
as some observed globular clusters show different popula- 
tions separated by sharp differences in helium abundances. 
The HDFs of the models illustrated in Fig. [4] are shown in 
Fig. [S] We now briefly discuss this figure: 

Panel a - This panel shows that in the RMMPRM the vast 
majority of the stars have Y ~ 0.25, with a thin tail extend- 
ing up to Y — 0.37. For smaller amounts of pristine gas 
(dashed histogram) the left peak is lowered and shifted to- 
ward larger values of Y (for the same reason, the cusp of the 
dashed line in panel a of Fig. [4] is shifted compared to the 
RMMPRM cusp). The tail in the range 0.29 < Y < 0.33, 
as well as the secondary peak at Y — 0.37, increase their 
relative importance. 

Panel b - The temporal shift of pristine gas accretion pro- 
duces its major effects on the secondary, helium rich peak. 
If the accretion is delayed (RMTACM, dotted line), a larger 
number of stars form with the pure, He-rich, AGB ejecta, 
giving rise to a non- negligible secondary peak at Y — 0.37. 
With an early accretion (dashed line) this peak is substan- 
tially erased because only few stars form with such an high 
Y value as the system quickly accretes pristine gas forming 
the main peak at Y ~ 0.26. Later on, as discussed in the 
previous subsection, most of the pristine gas is depleted by 
the SF and the ISM is composed of a larger fraction of He- 
rich AGB ejecta; thus the stars forming at this stage give 
rise to the HDF secondary peak at Y ~ 0.32. 

Panel c - From panel c of Fig. [4] we have seen that in the 
model with larger r (RMTAUM, dotted line) the cusp is 
slightly brought forward; for the same reason the main peak 
of the RMTAUM HDF is shifted a little "rightward" . The 
greater height of the secondary peak at Y — 0.37 in the RM- 
TAUL (dashed histogram) is due to the longer time during 
which the stars form from pure AGB ejecta, before the ac- 
cretion becomes substantial. 

Panel d - The RMEFFL (dashed histogram) creates little 
less O-poor, He-rich stars at early times, and the HDF shows 



therefore a slightly lower secondary peak at Y — 0.37 A 
smaller v leads to a longer persistence of the pristine gas 
and of its effect on the ISM chemical characteristics. For this 
reason the cusp of the dashed curve is shifted "upward" in 
panel d of Fig. [4j and the dashed HDF is shifted "rightward" 
in the present panel. 

Panel e - The models shown in panel e differ by the evo- 
lutionary time. Although the chemical paths are identical, 
the HDFs are quite different and reflect the different stellar 
crowding of the two arms (cf. panel e of Fig. |4|. The RMT- 
NDL histogram (shorter evolution) shows two well defined 
populations: a minor one with Y > 0.35, and the major one 
peaked around Y = 0.26. In the RMTNDM histogram a 
third population is present at Y = 0.31 as a consequence of 
the longer permanence of the model on the rising branch. 

Panel f - Given the larger amount of AGB ejecta in the 
RMFGDM (red dotted line), its He-rich peak is higher than 
in the RMFGDL. The dotted histogram also shows an in- 
termediate population in the range 0.3 < Y < 0.35 which 
is substantially absent in the dashed histogram. This inter- 
mediate population is a consequence of the bending toward 
lower values of O of the rising branch of the RMFGDM in 
the [O/Fe]- [Na/Fe] plane. 

5.2.3 Delayed star formation 



As pointed out in 12.11 in some clusters such as ui Cen 
and NGC 2808 multiple discrete populations are present, 
each characterized by a specific helium abundance. This led 
Rcnzini (2008) to argue that helium-rich material was accu- 
mulated in the ISM for a sufficient long time until suddenly 
a burst turned a major fraction of the ISM into stars. Ac- 
cording to this author, a continuous SF proceeding along 
with the ISM helium enrichment would have resulted in a 
continuous distribution of helium abundances in the newly 
formed stars, hence in a broadening of the MS rather than 
in well separated sequences as actually observed. 

In all the models we discussed so far the SF starts "ab 
initio" , as soon as the ejecta of the massive AGB stars begins 
to collect within the GC. To investigate the effect of a SF 
burst we run two models in which the gas cumulates in the 
GC since the beginning as usual, but with the SF inhibited 
up to the time £b,7 = 5 (RMTSFL, blue dashed line) and 
tb,7 = 7 (RMTSFM, red dotted line). The chemical patterns 
are very similar (see the lower panel in Fig. |S]), but with 
the dotted cusp shifted 0.1 dex "leftward" and "upward" 
compared to the dashed one. In fact, in the RMTSFM a 
large amount of AGB ejecta gathers before the occurrence 
of the SF; the dilution by the pristine gas is therefore less 
effective and the cusp takes place earlier. In this latter model 
the stars are essentially distributed along the rising branch, 
as expected given the large SF delay. 

In RMTSFM the effects of the pristine gas are less effec- 
tive, and the primary peak of the HDF (see the upper panel 
in Fig. [6]) is higher and shifted rightward (reflecting the shift 
of the cusp of dotted line in the lower panel of Fig. [6j. The 
secondary peak is absent because the SF was not active in 
the beginning, when the gas was very helium-rich. This sim- 
ple experiment shows that, with a SF delay large enough, it 
is possible to obtain a sharp intermediate population, but to 
the expense of the extreme one. Both populations can not 
be achieved with this mechanism. 
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Figure 5. Stellar helium distribution function for different models. The models in each panel are the same shown in the corresponding 
panel of Fig. [4] 



6 COMPARISON WITH REAL CLUSTERS 

In this section we apply the chemical evolution framework 
described in the previous section to build two specific mod- 
els for NGC 2808 and M4. These two clusters are indeed 
good prototypes to exemplify the cluster to cluster differ- 
ences. NGC 2808 is a complex mixture of populations, one 
of which is very extreme, and has a very high helium content 
testified by the presence of the blue main sequence. M4 is 
much less massive, and shows milder variations in elemental 
abundances, particularly in oxygen, but it apparently has a 
larger spread in sodium. It is to be preferred to other sam- 

! iles also because a ve ry accurate recent spectroscopic survey 
Marino et"al]|2008rj ) 



6.1 NGC2808 

The cluster NGC 2808 has a mass M ~ 1.6 x 10 6 M Q , 
and harb ours three main p o pulat ions, as initially specu- 
lated by iD'Antona &: Caloil i|2004h and successively con- 
firmed by IPiotto et all ( 20071 ). Each population has a dis- 
tinct helium abundance: ~ 50% of the stars are credited with 
a primordial helium abundance (Y = 0.25), ~ 32% with 
Y ~ 0.30 and ~ 18% with Y ~ 0.37 jCarretta et alj|2009d : 



ID'Antona fc Caloil 12008). No iron abundance difference is 
present in these populations, but they show an ex tended 
[0/Fe]-[Na/Fe] anticorrelation (|Carretta et al.ll2006T ). 

This anticorrelation is apparent in the bottom-left panel 
of Fig. [7] which includes data obtained with the ESO high- 
resolutio n multifiber spe c trograp h FLAMES/GIRAFFE at 
VLT by ICarretta et all (|2009cT l and data obtained with 
the hi gher resolution UVES spectrograph (|Carretta et all 
l2009bl ). 

We simulated the GC NGC 2808 with a model 
with the parameter set (iac,7,£<md,7,T7,p P r,p*,FG, v, x) = 
(5,6,0.3,2.36,241,1,0.5). The choice of x = 0.5 is done in 
order to get a current 50% of FG stars, as suggested by the 
helium abundance; this produces a more realistic compari- 
son with the observations. 

The time evolution of the ISM chemical composition 
is illustrated in Fig. [8] as shown in this figure, owing to 
the relative shortness of the evolutionary time the C+N+O 
abundance remains essentially constant in this model. We 
remind that the absence of a substantial CNO increment is 
an essential constrain for any model describing the chemistry 
of the majority of GCs (cf. sect. 12.4) ). 

The cyan circles in the bottom-left panel of Fig.[7]repre- 
sent a statistical sample of the SG stars obtained integrating 
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dashed line indicates the evolution of RMNOPG. 



the equation system in sect. [4] In populating the [O/Fe]- 
[Na/Fe] plane with such circles, we introduce a stochastic 
error in the range 0-0.1 dex, comparable to the observational 
errors, scattering the position of the stars. 

The blue circles represent the FG stars. Since the FG 
stellar population is assumed to be present "ab initio" in 
our model, these stars are not an outcome of the model, but 
are included in the figure to ease the comparison with the 
observational data. The pristine gas from which these stars 
formed is assumed to have [O/Fe]=0.4 and [Na/Fe]=0 and 
also in this case we introduce a stochastic error in the range 
0-0.1 dex in the abundances. The relative number of FG and 
SG stars is given by the value of th e param eter x. Following 
the terminology of ICarretta et al.l l)2008l ) , the blue circles 
reproduce the Primordial population, while the cyan ones 
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Figure 7. Results of our model for NGC2808. The following set of 
parameters has been assumed: (i a c,7, tend, 7, T, p pr ,P*, FGj v, x ) = 
(5, 6, 0.3, 2.36, 241, 1, 0.5). The top-left panel and the bottom-left 
panel show the [Mg/Fe]- [Al/Fe] and the [0/Fe]-[Na/Fe] diagrams 
relative to NGC 2808, respectively. The black squares and arrows 
are data obtained with the ESO high-resolution multifibcr spec- 
trograph FLAMES/GIRAFFE at VLT bv lCarretta et al.l J2006F ). 
Magenta triangles and arrows a re data obtained with the higher 
resolution UVES spectrograph l|Carretta et al.|[2009bl '). The blue 
and cyan circles represent a sampling of the FG and SG stars 
by our model, respectively. The red dashed line represents the 
gas trajectory within the diagram; the sampled stars in principle 
would be located on this line, but we introduced a random scat- 
ter in the range 0-0.1 dex in their coordinates in order to mimic 
the observational errors (see text for more details). The top-right 
panel illustrates the stellar helium distribution function. In the 
bottom-right panel the evolution of the following quantities is re- 
ported: total amount of gas (bold solid line), stellar ejecta (solid 
line), pristine gas (long-dashed line), SG stars (dashed line). 



give the Intermediate and Extreme ([0/Fe]< —0.4) popu- 
lations. We emphasize that without the modification in the 
chemical composition of the FG stars discussed above (cf. 
sect. 12. 2p the Extreme population and part of the Interme- 
diate ones could not be reproduced. In any case, stars with 
[Na/Fe]> 0.5 remain unexplained, and suggest that a fur- 
ther parameter is at work. In order to obtain a better fit of 
the data points, we could have increased the sodium pro- 
duced in the mass range 6.5 - 9 M© to [Na/Fe]~0.7. This 
procedure has two negative implications: first, it can not 
explain, in any case, the six points at [Na/Fe]>0.7 having 
upper limits for the oxygen values. In addition, such an ex- 
trapolation, implying an increase of the sodium yield with 
increasing initial mass, is not compatible with our present 
understanding of the super- AGB evolution (see Sect. [2~Tj) . 
We prefer then to postpone a better solution of this prob- 
lem to when more data (e.g., when the oxygen abundance 
for these stars will be measured) and better models will be 
available. The imperfect reproduction of these (few) data 
is not an argument to believe that the origin of the gas 
forming these SG stars must be looked for somewhere else 
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Figure 8. Time evolution of the chemical abundances in the ISM 
of the NGC 2808 model. 



Figure 9. Time evolution of the chemical abundances in the ISM 
of the M4 model. 



than among super-AGBs. In fact, no other source can be 
able to deal, at the sam e time, with the repro duction of the 
Mg-poor, Al-rich stars (JDecressin et al]|2007l ). and with the 
extreme but homogeneous helium cont ent that is necessary 
to reproduce the triple main sequence (|D'Ercole et al.ll2008l ; 
iD'Antona fc CaloJl20oilRenzinill2008l . and Fig.0 top right 
panel). In summary, we regard the fit of NGC 2808 data as 
very satisfactory, in spite of this remaining problem. 

The three [Al/Fe]-rich stars in the upper-left panel 
of Fig. [7] coincide with the three magenta arrows in the 
[O/Fe]- [Na/Fe] diagram . In both diagrams these UVES data 
l|Carretta et al.l l2009bh show a very "clumpy" anticorrela- 
tion, with only Primordial and Extreme populations present, 
and a gap in the middle where no star is found. On the other 
hand, our model, giving rise to stars occurring throughout 
the whole O-Fe anticorrelation span, necessarily produces a 
"continuous" distribution of stars in the [Mg/Fe]-[A1/Fe] di- 
agram. We believe that the gap in the Mg-Al anticorrelation 
displayed by the UVES data is not real (as indeed sh own by 
the FLAMES/GIRAFFE data (ICarretta et al.ll2009ch in the 
[O/Fe]- [Na/Fe] plane), and we predict that an Intermediate 
population in the [Mg/Fe]-[A1/Fe] diagram will be observed 
in the next future. 

Finally, the top-right panel of Fig. shows the helium 
stellar distribution. Three stellar populations are clearly 
present. The most He-rich stellar group (Y > 0.35) origi- 
nates from the FG ejecta before the pristine gas enters the 
GC, while the Intermediate population (0.3 < Y < 0.35) 
forms after the accretion of the original gas. The ratio be- 
tween Extreme and Intermediate He-rich stars is 0.56, as 
observed. Contrary to the abundance of the FG stars, which 
is simply imposed, the value of the ratio between Extreme 
and Intermediate stars is a genuine result of our model. 



6.2 M4 (NGC 6121) 

The globular cluster M 4 has a mass M = 6.3 x 10 4 M 
l|Mandushev et alJ Il991f ) and its color-magnitude diagram 
does not show any of the fingerprints of the presence of 
multiple populations found in the most massive clusters. 
However, the O-Na anticorrelation is present in this clus- 
ter as in all the Galactic GCs. Instead, cont rary to NGC 
2808, t here is no Mg-Al anticorrelation in M4 (|Marino et al.l 
l2008bT ). As the stellar evolution must be the same in all the 
GCs, our aim is to fit M4 data adopting the same chemical 
composition of the AGB ejecta utilized in modelling NGC 
2808. However, we must assume a different composition of 
the pristine gas. This can be understood by looking at the 
observational data in Fig.lTland llOl The FG oxygen, sodium, 
magnesium and aluminium for NGC 2808 are [0/Fe]~ +0.4, 
[Na/Fe]~ 0.0, [Mg/Fe]~ +0.28 and [Al/Fe]~ +0.15. In M4 
we have [0/Fe]~ +0.5, [Na/Fe]~ +0.1, [Mg/Fe]~ +0.5 and 
[Al/Fe]~ +0.4. Therefore, with respect to the input values 
for NGC 2808, in the present model we adopt O and Na 
abundances 25% higher, and Mg and Al abundances 67% 
and 78% higher, respectively. As a consequence, we should 
also consider that the initial abundance of 20 Ne, another a 
element, may have been larger in the pristine gas. Follow- 
ing the test computations described in section 2.2, also the 
sodium yield of the AGBs will be larger, and we can rea- 
sonably increase by a factor 0.2 dex the ratio [Na/Fe] of the 
ejecta of the AGB of any mass. Notice that this last change 
in the yield table is only applied to the M4 case, due to the 
different initial chemistry of its FG. 

For this model of M4 we set the pa- 
rameters (tac,7,i C nd,7,T7,Ppr,P*,FG,l / , x) = 

(6.5,10.5,0.91,9.4,0.1,0.7). The results are illustrated 
in Fig. [9] and Fig. 1101 The former figure shows the evolution 
of the ISM chemical abundance of several elements; in 
particular, it can be seen that the variation of C+N+O 
abundance is not larger than 0.1 dex, as required by the 
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Figure 10. As Fig. [7] but for M4. Data are from iMarino et al.l 
(|2008bh . In this case, for a better comparison with the data, the 
error on [Na/Fc] is taken 0.05 dex instead of 0.1 dex as in the 
case of NGC 2808. 



observations IjYong et al.ll2009r ). Figure [10] can be compared 
to the analogous Fig. [7] relative to NGC 2808 to understand 
the different role played by various parameters in the two 
cases. 

Our model shows that the relative amount of pristine 
gas must be larger than in NGC 2808. The reason is illus- 
trated in the bottom-left panel of Fig. 1101 The Na-O anti- 
correlation is less extended than in NGC 2808; in particular, 
all the stars in M4 exhibit a positive ratio [O/Fe], while in 
NGC 2808 stars are found with ratios as low as [0/Fe]~ —1. 
In order to reduce the fraction of this kind of stars we set 
v = 0.1 so to minimize the SFR in the beginning, before the 
accretion of pristine gas, when the ISM is composed only of 
O-poor AGB ejecta. After the pristine gas accretion and the 
consequent increase of the ISM density, the SFR becomes 
significant despite the low value of v. Thus, a large fraction 
of SG stars form from an ISM strongly "diluted" by the 
pristine O-rich gas. 

In order to form the observed Na-rich stars, the model 
for M4 also requires a value of £ on d larger than that of NGC 
2808. By extending the SF duration, Na-rich stars can form 
at late time when the pristine gas has been strongly depleted 
by SF and the ISM is mainly composed of the Na-rich stellar 
ejecta of stars with M < 6 Mq . 

The chemical evolution of the GC is summarized by 
the dashed line in the bottom-left panel of Fig. 1101 The 
model starts at the O-poor extreme and moves rightward 
relatively fast and in a low SFR regime, so that only ~ 3% 
of the total SG population form during this excursion. The 
trajectory has a cusp at [0/Fe]~ 0.4 where then develops a 
rather vertical branch due to the reduction of pristine gas 
and the increment of Na-rich ejecta in the ISM "mixture" . 
The longer length of time spent along this branch and, more 
importantly, the larger SFR during this stage determine the 



clustering of the SG stellar population in this area of the 
Na-O plane. 

The amount of pristine gas accreted by the GC is a key 
parameter to fit the observed data both in the [0/Fe]-[Na/Fe] 
and the [Mg/Fe] -[Al/Fe] diagrams. An increment of pristine 
gas causes the abundances of SG stars to approach those 
of FG stars and, in particular, reduces the differences in the 
[Na/Fe] and [Al/Fe] ratios between Primordial and Interme- 
diate stars. The Na abundance is, however, more sensitive 
to this effect: a model of ours (not shown here) reveals that, 
after trebling the accreted gas, the Al abundances are still 
consistent with the data, while stars with [Na/Fe] > 0.4 can 
not be reproduced. 

Finally, from the top- right panel of Fig. [I0]we note that 
the helium abundance of the SG Intermediate population is 
not much larger than that of the FG one, and the two popu- 
lations, although still distinguishable, blend together in the 
HDF. This small helium spread can be revealed in care- 
ful observat ions of the main se quence width (see, e.g., the 
analysis by iMilone et al.ll2010l . for the cluster NGC 6752) 
and for its possible conseque nces on th e horizontal branch 
(HB) morphologies (JD'Antona fc Caloil l200ct) . The HB of 
M4 shows signs of bimodality (jNorria Il98lh but a careful 
representation of the data is very tricky, as it requires use 
of models including also the CNO abundance enhance ments 
of the SG (see the discussion by ISalaris et al.l [2008, con- 
cerning the bimodal d istribution in the HB of NGC 1851). 
iD'Antona et al.l (120091 ) notice that the blue side of the HB of 
M4 might be slightly overluminous with respect to the blue 
HB of NGC 1851, and tentatively attribute this feature to 
possible helium enhancement in these stars. This would be 
consistent with our helium distribution in M4. 

The Extreme population is still present, but is only 2% 
of the whole GC stellar population. As observationally this 
population is lacking in the clusters that are not very mas- 
sive, this may be an indication that the formation of the 
Extreme population is inhibited below some critical mass. 
Alternatively, given its paucity, we can guess that this pop- 
ulation exists but is not yet observeqj. 

In any case, as shown in sect. 15.21 it is enough to change 
somewhat one or more parameters to create or remove small 
structures in the helium distribution function. Thus, at this 
stage, no firm conclusion can be drawn on the Extreme pop- 
ulation we found, or on the detailed helium spread of the 
Intermediate population. More realistic models are planned, 
in which the chemical machinery developed in this paper 
will be implemented in a hydro-code. 



7 CONCLUSIONS 

In this paper we have studied the origin of the chemical 
patterns which have been observed in many globular clusters 
and which are considered the spectroscopic fingerprints of 
the presence of multiple stellar populations. Specifically, in 
our investigation we have focussed our attention on the O- 
Na and Mg-Al anticorrelations and the helium distribution 
function. 



6 Notice however that a very helium rich, even small, population 
would be photometrically very evident, as it would popu late the 
extreme blue horizontal branch llD'Antona fc C aloi 2008). 
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In our model we have assumed that AGB stars are the 
polluting source of gas with the anomalous abundances of 
light elements observed in globular cluster second-generation 
stars. Our chemical framework is based on a one-zone model 
following the formation of SG stars from a mix of ejecta of 
AGB stars and pristine gas. 

We have carried out a large number of simulations to 
explore the dependence of the SG chemical properties on the 
parameters characterizing the star formation process and 
the dynamics of the involved pristine gas. Finally, we have 
used our framework to model the observed chemical patterns 
in the Galactic clusters NGC 2808 and M4. 

The main results of our study are the following: 

(i) The current stellar models provide AGB ejecta in 
which both Na and O decrease with increasing stellar mass. 
In order to reproduce the observed O-Na anticorrelation, the 
gas from which SG stars form must be diluted with O-rich, 
Na-poor pristine gas (see e.g. Figs. [2] and [3]). 

(ii) The amount of pristine gas involved in the SG for- 
mation process, the timescales driving the dynamics of such 
gas, and the star formation efficiency play a key role in de- 
termining the extension of the O-Na anticorrelation and the 
fraction of extreme Na-rich/O-poor stars. 

(iii) The helium abundance distribution function is cor- 
related with the distribution of stars in the O-Na plane. Ex- 
treme Na-rich/O-poor stars are also those with an extreme 
He enrichement. Our models predict that all the clusters 
with a very extended O-Na anticorrelation should also host 
a population of He-rich stars. 

(iv) Our models show that the extension of the O-Na an- 
ticorrelation is closely correlated with that in the Mg-Al 
plane; Na-rich/O-poor stars have also high Al and low Mg 
abundances. 

(v) We have used our framework to build specific models 
for two prototypical Galactic globular clusters: NGC 2808, 
a massive cluster that hosts a SG population characterized 
by a very extended O-Na anticorrelation and includes a very 
He-rich population, and M4, a low-mass cluster with a sig- 
nificantly less extended O-Na anticorrelation that does not 
include extremely Na-rich O-poor stars and for which there 
is no photometric evidence of a He-rich population. Despite 
the significant differences in their chemical patterns, in both 
clusters a significant f raction of stars belon g to the SG (50 
per ce nt in NGC 2808 , [Ca"rretta et al. 2009c, and 65 per cent 
in M4 lCarrettaet all |2009ci ): [Marino et all (|2008bl ))'). 

Our models successfully reproduce the differences in the 
O-Na anticorrelation observed in these two clusters, the dis- 
tribution of stars in the Mg-Al plane in M4 and predict an 
extended Mg-Al anticorrelation for NGC 2808. Al and Mg 
abundances for NGC 2808 have been determined with UVES 
observations only for a small number of either Na-rich/O- 
poor or Na-poor/O-rich stars. As predicted by our model, 
these stars populate only the Al-rich/Mg-poor and the Al- 
poor/Mg-rich regions of the Mg-Al plane. According to our 
model, future determinations of Al and Mg abundances of 
targeting stars with intermediate Na and O should lead to 
intermediate Al and Mg abundances filling the extended Al- 
Mg anticorrelation. 

(vi) The helium distribution, including the extreme pop- 
ulation formed directly from super-AGB ejecta is well repro- 
duced by our model for NGC 2808 (see Fig. [7|. In clusters 



like M4, in which a larger dilution with pristine matter is 
necessary to model the O-Na and Mg-Al patterns, a large 
helium dispersion is not required, but a small helium spread 
should generally be present (see Fig. I10[) . 

Our investigation has shed light on the key chemical 
and hydrodynamical ingredients determining the formation 
of the chemical patterns observed in globular clusters. The 
results presented here are to be the starting point inform- 
ing further study based on full hydrodynamical simulations. 
Several issues will require further investigation. Specifically, 
as for the stellar evolution models, we have shown that our 
conservative choice of an educated linear and monotonic ex- 
trapolation for the sodium and oxygen abundances in the 
mass range from 6.5 to 9 Mq reproduces the general trend 
of the O-Na ant icorrelat ions of different clusters, when also 
the differences in neon abundances in the FG of different 
clusters are taken into account. In NGC 2808, however, we 
can not reproduce the (few) very large sodium values of some 
stars for which only upper limits on the oxygen abundance 
are available (see Fig. [7J). Another parameter is probably 
at work here, and only further investigation into the super- 
AGB phase will shed light on this issue. 

The lack of any evidence of a significant metal enrich- 
ment in most clusters hosting multiple stellar populations 
implies that neither ejecta from FG nor SG supernovae is 
involved in the chemical enrichment of matter from which 
SG stars form. While in our model SG formation starts after 
the end of the FG SN II epoch, it is to be further clarified 
whether the lack of metal enrichment from SG supernovae 
is due to SG forming with IMF truncated at M < 9 M 
or whether this is a consequence of a more complex gas 
dynamics to be further explored with full hydrodynamical 
simulations. 

Finally, additional full hydrodynamical simulations 
will be needed to clarify the source of pristine gas 
and the accretion mechanism. Different processes have 
been considered in the literature, su ch as accretion from 
a diffuse surrounding medium (e.g . iLin fc Murray! 120071 : 



IPflamm-Altenburg fe Kroupa 20091) or interaction with 



molecular clouds ( Bekki fc Mackevll2009r ) . Conditions for an 



effective gas collection within the cluster turn out to be 
rather specific. Instead, our model requires a mechanism 
general enough to work in all the clusters and able to take 
into account the cluster-to-cluster differences in the dynam- 
ics and the amount of pristine gas involved in the SG for- 
mation process. 
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